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ABSTRACT 

The gravitational evolution of the cosmic one-point probability distribution function (PDF) has been 
estimated using an analytic approximation that combines gravitational perturbation theory with the 
Edgeworth expansion around a Gaussian PDF. Despite the remarkable success of the Edgeworth ex- 
pansion in modeling the weakly non-linear growth of fluctuations around the peak of the cosmic PDF, 
it fails to reproduce the expected behaviour in the tails of the distribution. Besides, this expansion is 
ill-defined as it predicts negative densities and negative probabilities for the cosmic fields. This is a 
natural consequence of using an expansion around the Gaussian distribution, which is not rigorously 
well-defined when describing a positive variate, such as the density field. Here we present an alterna- 
tive to the Edgeworth series based on an expansion around the Gamma PDF. The Gamma expansion 
is designed to converge when the PDF exhibits exponential tails, which are predicted by Perturbation 
Theory, in the weakly non-linear regime, and are found in numerical simulations from Gaussian initial 
conditions. The proposed expansion is better suited for describing a real PDF as it always yields positive 
densities and the PDF is effectively positive-definite. We compare the performance of the Edgeworth 
and the Gamma expansions for a wide dynamical range making use of cosmological N-body simulations 
and assess their range of validity. In general, the Gamma expansion provides an interesting and simple 
alternative to the Edgeworth series and it should be useful for modeling non-gaussian PDFs in other 
contexts, such as in the cosmic microwave background. 



1. INTRODUCTION 

We aim at studying one dimensional probability density 
functions (PDF), p(5), which characterize the statistical 
properties of a stochastic field at a single point 5 = S(r). 
Here we shall concentrate on fluctuations of the density 
field: 5 = p/p — 1 (being p the mean value of the den- 
sity field, p), smoothed over some fixed scale R. However 
many of the arguments presented here are quite generic 
and could also be applied to other contexts. The approach 
we will follow is to try to recover the full shape of the PDF 
from the knowledge of its first order moments. This has 
become by now a textbook problem and there are several 
standard ways to address it (for a review see Kendall, Stu- 
art & Ord 1994). The solution is not unique, unless there 
is a well defined hierarchy in the moments and we can 
define some perturbative approach to the problem. 

There are a number of studies to predict the evolution 
of clustering of density fluctuations, and in particular of 
its PDF. There have been attempts to derive the PDF 
from analytic approximations, such as the Zeldovich Ap- 
proximation (Kofman et al. 1994). Although the Zel- 
dovich Approximation reproduces important aspects of the 
non-linear dynamics, it is a poor approximation for the 
PDF and its moments. One way to improve that is to 
take advantage of the exact non-linear perturbation the- 
ory (PT) to estimate the moments (Bernardeau 1992) and 
use them to derive the PDF from the Edgeworth expansion 
(Juszkiewicz et al 1995, Bernardeau & Koffman 1995). In 
this case the PDF is predicted to an accuracy given by the 
order of the cumulants involved. The Edgeworth expan- 



sion has since been used as a tool to characterize the PDF 
of matter (eg Kim & Strauss 1998, Blinnikov & Moessnew 
1998, Pen 1998) and CMB fluctuations (Amendola 1996, 
Popa 1998). Earlier phenomenological approaches to the 
construction of the cosmic PDF were developed by Saslaw 
& Hamilton (1984) and Coles & Jones (1991). 

One serious shortcoming of the Edgeworth approach is 
that the series yields a PDF that is ill-defined. It has neg- 
ative probability values and assignes non-zero probability 
to negative densities (5 < —1). This latter problem origi- 
nates on the fact that the Gaussian PDF, which is the basis 
for the Edgeworth series, only makes physical sense when 
the rms fluctuation a is very small. Here we will try to ad- 
dress some of these problems by exploring the possibility 
of carrying out expansions around better behaved PDFs, 
more suitable to yield positive densities when the variance 
is not that small. We shall concentrate on the Gamma 
PDF, but other distributions may be handled within the 
same framework, as our general analysis will show. 

Whenever the moment generating function of the PDF is 
known, it is then possible to reconstruct the full PDF. This 
has been done previously by using the Legendre trans- 
form (Fry 1985) or the inverse Laplace transform (Balian 
& Schaeffer 1989, Bernardeau 1992, Bernardeau & Kof- 
man 1995). Since for gravitational clustering in the weakly 
non- linear regime the variance of the PDF is small, one can 
expand the above transforms to recover the PT limit. In 
the case of Gaussian initial fluctuations this perturbative 
expansion yields to the well-known Edgeworth series. 

In §D we shall introduce the Edgeworth series as a 
generic saddle point approximation to the PDF, along the 



lines of Fry (1985). An alternative expansion, around the 
Gamma distribution, is introduced in §||. It is expressed 
as an orthogonal polynomial expansion of an arbitrary 
PDF in terms of Laguerre polinomials. The latter are 
the counterparts to the Hermite polinomials when one ex- 
pands around an exponential tail instead of a Gaussian 
one. A detailed comparative analysis of the Edgeworth 
and Gamma expansions with respect to N-body simula- 
tions is presented in §|] and §^[ A final discussion with our 
conclusions is given in §p[ 

2. EXPANSIONS AROUND A GIVEN PDF 

2.1. One-Point Statistics 

As usual, we shall denote statistical averaging by brack- 
ets: (•■■), so that the expectation value for the moments 
are: 



mj = (5 J )= / P (8)5 J d6 



(1) 



with J an integer that labels the order of the correspond- 
ing moment. J = 1 corresponds to the mean, which for the 
density fluctuation is zero, mi = ( S ) = 0. In this notation 
the variance, Var (<5), and rms fluctuation a, are defined 
as: Var (S) = a 2 = m 2 — rn\. It is useful to introduce the 
cumulants kj: 



(6 J ) C 



d J ^(t) 



dt J 



d J log M(t) 



t=o 



dt J 



(2) 



t=Q 



where ip(t) = \ogA4(t) is given in terms of the moments 
of the PDF through M(t): 
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Gravitational clustering from Gaussian initial conditions 
predicts < 5 J > c ex < 5 2 > c on large-scales, thus it is 
more convenient to introduce the following ratios, 
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where the Skewness, S3, is the third-order ratio, and the 
Kurtosis, Si, is the fourth-order one. From ip(t) we can 
get back the PDF, p(S), by using the inversion formula: 



p(5) = 



-\-ioo 



dt 

27 
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(5) 



Consider two differentiable PDFs, pi{5) and P2(5), with 

cumulants kj and kj , it follows that (see eg. Kendall, 
Stuart & Ord 1994): ' 



p^'(5) = exp 
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This equation is easy to prove by reobtaining the moments 
of pW through the generating function, having assumed 



that those of p^ are given by k j . After partial integra- 
tion, the generating function of p^ 1 ' arises immediately, 
what proves the above equality (it is a nice exercise). 

Equation (g) allows one to use the most convenient PDF 
to do the series expansion, p^ 2 >{8). In particular, if one 
uses the Gaussian as the parent PDF one ends up with 
the Gram-Charlier series which yields the Edgeworth ex- 
pansion in the perturbative limit (ie, when the variance of 
the PDF is small). 

In this paper, we will focus on the Gamma distribu- 
tion. In that case, we have that p( 2 '(5) is given by the 
Gamma PDF (see Eqdjl) below) with S = (3z + a. Naming 
pW(<5) =pW( z ) (i = 1,2), we get 
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However such a straightforward approach is ill-defined as 
the measure of the expansion is proportional to the order 
considered (the order of the expansion, J, is involved in 
the index denotingthe order of the generalised Laguerre 
polinomials in Eq.[Q). 

This will force us to modify the general method to derive 
a consistent expansion in terms of the relevant orthogonal 
polynomials (see §0). 

2.2. The Gaussian vs. the Gamma PDF 
In the case of the Gaussian (or normal ^(0, er)) PDF: 



p{5) = p G {5) = 



1 



V2 



exp 
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where a is the rms standard deviation, the only parameter 
in this distribution. As the overdensity has to be positive 
p > 0, we have that S < — 1, and a Gaussian PDF only 
makes physical sense when a — > 0. For the Gaussian, 
Sj = kj = for J > 2. 

The Gamma — also called negative binomial or Pear- 
son Type III (PT3) PDF— arises from the Chi-Square 
distribution with N degrees of freedom when 1/er 2 = N/2 
is taken to be a continuous parameter. This yields 



P (S) = 4>{S) 
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(9) 



The hierarchical coefficients in Eq(|]) are constant for all 
values of the variance in this case and give Sj — (J — 1)!. 
These Sj values are equal to those of a simple exponential 
distribution. The Gamma PDF (or similar ones) has been 
found to be useful at describing the galaxy distribution 
(see eg, Fry 1986, Elizalde & Gaztahaga 1992, Gaztahaga 
& Yokohama 1993). 

2.3. The Saddle Point Approximation 

The moment generating function summarizes all the in- 
formation concerning the higher-order cumulants, as long 
as the series expansion in terms of the latter converges. In 
the majority of the cases this is true (one counterexam- 
ple to this rule is the Lognormal distribution). Thus, one 
may reconstruct the PDF from the the moment generating 
function in a consistent way as shown in Eq(pf). 



To obtain an asymptotic expansion of p(S) for small 5, 
we introduce the Legendre transform, 



S = dip{t)/dt, G(d) = 8t-ip(t) ) 



(10) 



where the convexity of G(5) is related to that of ip(t). Re- 
placing this in the original expression for p(S), we are left 
with 



p(S) 



G'=+ioo G ,, d5 
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exp[SG'(5) + SG'{S) - G(S)] 



(11) 

which is_dominated by stationary points of the exponential 
at 6 = S (for real finite (5). The Saddle Point approxima- 
tion of this integral is given in Fry (1985) and follows from 
the usual approach (Morse and Feshbach 1953): 



p(6) ~ [G"(6)/2*]V 2 exp[-G(£)]. 
Once normalized, the distribution reads as 
[G"(6)/2tt}^ exp[-G(<5)] 



p(S) 



fI™[G"(S)/2*] 1/a e*P[-G(S)]M 



(12) 



(13) 



Provided one constructs the generating function i/j(t) out 
of the irreducible moments, 

oo 1 1 1 

m = ££*- = \^ 2 + l s ^ 3 + k Si(jHi + 

n= 2 

(14) 
(where the latter equality shows the expansion in terms 
of the hierarchical amplitudes) , there is a general develop- 
ment for G(6) in powers of S which is derived by inverting 
the t variable in the Legendre transformation as to get 
t = t (5), leading to 



G(6) 



-6 2 ^S 3 
2 6 



- f S 3 2 ~ § I * 
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(15) 

To get a proper expansion around the Gaussian PDF, we 
first need to arrange the exponential by factoring out the 
quadratic term. In the limit of small a, we then get 



exp[-G(i/)] 



1 + ±S 3 v 3 a + I (5 4 
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(16) 



where v = 5/a. Applying the same expansion to the de- 
nominator (the normalization) and recalling the property 
for Gaussian integrals 
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we obtain a polynomial in powers of the small parameter 
a 2 of the form: 
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Finally, we multiply both developments derived above and 
keep terms up to a 3 , to end up with 



p(v) 
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H n (y) being the Hcrmitc polynomials, 
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which is the well-known (perturbative) Edgeworth series 
of a PDF up to third order. Higher-orders in the Edge- 
worth series can be obtained by keeping higher-orders in 
the Taylor expansions of Eqs(|l4|),(|l5). 

We stress that the latter expansion is derived under the 
assumption that the distribution is hierarchical, ie, the Sj 
are independent of a. According to this, the Edgeworth 
expansion may be generalized to non-hierarchical PDFs 
whenever the scaling of Sj(a) is known and replaced in 
the generating function ip(t) from which the Saddle-Point 
approximation of p(S) is built. 

Non-linear PT for Gaussian initial conditions predicts 
corrective (c dependent) terms to the leading order contri- 



bution to Sj of the form Sj 
? (o) c(l) 



S 



(o) 



■s<V- 



-0(er 4 ), where 



S\ ,Sj are coefficients independent of a. This must be 
taken into account to make consistent predictions from 
non-linear dynamics for the third-order term (or higher) 
in the Edgeworth series of the PDF (see Bernardeau & 
Kofman 1995). 

The Edgeworth expansion — or any other expansion 
based on the symmetry around the peak of the PDF be- 
ing approximated — is only a good candidate for fitting 
the evolutionary picture of the density profile as a first 
order approach, because high-density (non-Gaussian) ex- 
ponential tails develop in further stages of the non-linear 
evolution for arbitrary initial conditions. 

The Edgeworth series has also been applied to non-linear 
transformations of the Gaussian process to fit the exponen- 
tial tails observed in the simulations as the system evolves. 
The 'skewed' Lognormal approximation put forward by 
Colombi (1994) is an example of this scheme which takes 
advantage of the apparently privileged role the Lognor- 
mal PDF plays among the non-Gaussian ones. This is 
suggested by the integration of the continuity equation in 
Lagrangian coordinates on one hand, and by the good fit 
to the observational PDF on intermediate scales (related 
to an — 1 spectral index for a power-law power spectrum) , 
on the other hand. 

Nevertheless, it is still lacking in the literature a well- 
defined expansion around a non-Gaussian PDF which may 



be better suited than the Gaussian to model the gravita- 
tional evolution of cosmic fluctuations in the weakly non- 
linear regime. It is not clear yet whether the initial condi- 
tions of structure formation in the universe were Gaussian 
(as suggested by standard inflationary models), or not. In 
the latter case, it is necessary to investigate expansions 
around non-Gaussian PDFs if one wants to describe clus- 
tering. These issues have been our main motivation to in- 
troduce an alternative expansion around a non-Gaussian 
PDF. 

In the next section we derive a general expansion around 
the Gamma distribution, ie., around an arbitrary exponen- 
tial tail, making use of the completeness and orthogonality 
properties of the Laguerrre polynomials. They are for- 
mally analogous to the Hermite polynomials that appear 
in the Edgeworth series around a Gaussian. 

3. EXPANSION AROUND THE GAMMA PDF 

Our starting point here will be the expansion of the PDF 
in terms of the Gamma distribution, with a basis that will 
be given by generalized Laguerre polynomials (instead of 
the Hermite ones in Eq(piO|)). That such an expansion 
makes sense becomes clear from the fact that the Gamma 
distribution is proportional to the measure associated with 
this particular family of orthogonal polynomials. In this 
sense, the Gamma expansion is formally reminiscent of the 
Edgeworth series as the main difference consists of replac- 
ing the Hermite polynomials by the Laguerre ones as the 
basis for the expansion. Notice however that while the 
Edgeworth series might be built from the Gram-Charlicr 
series (Cramer 1946), given by the Gaussian PDF and its 
derivatives (see Juszkiewicz et al. 1995), here we cannot 
simply take successive derivatives of the Gamma distribu- 
tion. This is because, in doing so, we would actually fail to 
make contact with a consistent theory of orthogonal poly- 
nomials. In short, if the generalized Laguerre ones are to 
be used (and those are the only possibility in the case con- 
sidered), it turns out that the integration measure changes 
with the order of the generalized polynomial (and not just 
with the parameter of the family), and this would inval- 
idate the whole approach. We should, by the way, recall 
the good properties of a well defined expansion in terms of 
an orthogonal basis of polynomials, that arc orthonormal- 
ized with respect to an scalar product defined by a fixed 
integration measure. This is a rigorously defined mathe- 
matical theory, that in the case of the Gaussian measure 
happens to coincide with the Taylor expansion in terms of 
the derivatives of the function. In the case of the Gamma 
function, on the contrary, it turns out that the two ex- 
pansions do not coincide, and only the one in terms of 
orthogonal polynomials has rigorous mathematical sense. 
Thus, this must be the starting point in our approach, 
as the Taylor expansion (given in terms of the successive 
derivatives) must be relegated to a mere formal expansion 
lacking adequate convergence properties. 

The key point in our scheme is to build a general and 
consistent expansion in terms of the relevant orthogonal 
polynomials (the Laguerre ones, in the present case), not 
necessarily given by derivatives of the parent PDF (ie, we 
do not make use of Equation (0)). Bearing this in mind, 
we modify the approach based on the general Eq(ph and 
define an expansion of the PDF in terms of the Gamma 



distribution as follows: 

oo 



(20) 



ic, we define an expansion for which all orthogonal polino- 
mials are well defined, with coefficients: 

c ^rfr^r^ )4P " )(M) ^ (21) 

being 4>(fi) the Gamma PDF: 
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x — a 



(22) 
(23) 



This is actually a three-parameter (p, a and /3) family of 
distributions out of which only one, p, is relevant for nor- 
malized variables (such as density fluctuations, 5). 

The generalized Laguerre polynomials we shall need are 
given by: 



4 P - 1) (M)=E^ L ( n n + 4"> 
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in particular, 
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and the coefficients c„ are easy to calculate: 

c = 1, ci = c 2 = 0, 

r(p + i) 
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T(p + 5) 
In general, one finds: 
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It is clear that for Sk = (k — 1)! we recover the Gamma 
PDF, p{n) = <KM)- 

Note that for the first coefficients it holds 

fln,n-i = n(n- 1), . . . 



We should stress the fact that for higher orders (ie., for 
Cq or higher) there appear linear terms in p = 1/ct 2 which 
couple to those quadratic in [Sk — (k — 1)!]. For instance, 



<•'. = ^ + l.l [(Se - 5!) - 30(^5 - 4!) + 300(S 4 - 3!) 
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1200(S 3 -2!) + 10p(S 3 -2!) 2 ] 
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Notice that an appropriate expansion for small a will have 
contributions from these higher-order terms. In order to 
see clearly which is the parameter in the expansion (Effl, 
and set up a comparison with the Edgeworth expansion, 
we express Eq(|2G) in terms of the same variable, v = S/a, 
so that y, = l/cP + u/cr. Doing this, we have: 

^ n— 3 v ' ' 

x V(-l) fc a„,fe[^-(fc -1)!]1 <f>(y) 



fc=3 



= <1 



ra=3 



x Y(-l) n - k a Thk {Sk ~(k- 1)!] <j>{v), (29) 
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with the F n = F n (v, a) being of the form: 



F n (p, a) = -H n {v) + QJv) a + K n (v) a 2 + 0(a 3 ) (30) 
n! 

with H n the Hermite polynomial of order n, and 
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To summarize, we see that what we have in Eq(E9f) is in 
fact an expansion in terms of the Gamma PDF in the per- 
turbative limit (ie., when a is small), in formal analogy 
to the Edgeworth expansion (which uses the Gaussian as 
the parent PDF). The new expansion should presumably 
be much better suited than the Edgeworth expansion, to 
parametrize PDFs with exponential tails. 

4. COMPARISON OF THE EDGEWORTH & GAMMA 
EXPANSIONS 

We can now compare the Gamma with the Edgeworth 
series. We just have to replace Eqs(|30|),(pi|) in Eq(p9|), 
and express the (third order) Gamma expansion as a se- 
ries in cr, the rms standard deviation. We can write it in 
the following compact notation: 



PJ V ) 
4>{v) 



= 1 + H 3 (is)A 3 a 



+ {Hi{u)(A 4 - 3A 3 ) + H 6 (iv)A 2 3 /2 + 6 Q 3 (^)A 3 } a 2 
+ {H 5 {v)(A 5 - 4A 4 + 6A 3 ) + H 7 H(A 4 A 3 - 3A 2 ) 
+ H 9 (v)Al/6 + 24Q 4 H(-3A 3 + A 4 ) 
+ 360Q6MA 3 + 6^ 3 (^)A 3 }a 3 + O((7 4 ) (32) 

where <j){v) is the Gamma distribution (see Eq(g)) for the 
dimensionless variable v. The polynomials Q n ,7?.„ are 
given above, Eq(|3l]), and the reduced moments of the p{v) 

PDF, Sj, appear as differences with the moments, Sf of 
the parent PDF (over which we are expanding): 



A, 



S j — S 
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(33) 
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In particular, for the Gamma expansion, Sj = (J — 1)!. 
Eq. [gammaexp] is the main result of this paper. 

For the Edgeworth series around the Gaussian one has 
up to the same order, 
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where pg^) is the Gaussian PDF with Sf 

By comparing the Gamma and the Edgeworth expan- 
sions given above, we see that the Gamma expansion re- 
covers all the terms that appear in the Edgeworth expan- 
sion plus some corrective terms. In general the Gamma 
expansion has, by construction, exponential tails and a 
better general behaviour than the Edgeworth expansion, 
both on the positivity of p{y) and the variate itself, p. 

Notice that, as suggested by expansions Eqs(|32;),(|35|), 
the convergence of these series depends on the magnitude 
of the coefficients A,/ a J ~ 2 that weight the contribution 
from every polinomial. This is also true for the Edge- 
worth expansion, where the natural smallnes parameter is 
given by Sj a J ~ 2 (see also Juszkiewicz et al. (1995)). This 
can already be guessed from the Saddle Point approach to 
the Edgeworth series. In particular (see Eq(|l4[)), the cu- 
mulant generating function, ip(t), is expanded in terms 
of powers of at, with coefficients given by the cumulants 
kj = Sj (7 J ~ 2 . Accordingly, in the Gamma expansion, the 
convergence depends on how close the tail of the PDF is 
to an exponential one within the perturbative limit. 

Figure [I] (left panel) shows a comparison of the two 
expansions to leading order in cr, for a — 0.2 — 0.3 and 
S3 = 3 — 4, as labeled in the Figures. For reference we also 
show the Gaussian distribution with the same a (dotted 
line). The Edgeworth expansion (continuous lines) quickly 
develops negative probabilities (which are shown in abso- 
lute value) for negative values of v = 8 /a. It is also clear 
how the Edgeworth has Gaussian tails dropping quickly 
to zero, while the Gamma expansion (shown as a dashed 
line) has exponential-type tails. 

Figure p] also shows a comparison of the two expansions 
to the next order (second order) in a, for a = 0.2 — 0.3, 
S 3 — 4, and S4 = 15 — 20 as labeled in the Figures. 
The Edgeworth expansion (continuous lines) shows now 
more negative values of the probability (shown in absolute 
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Fig. 1. — Comparison of the leading order (left panel) and second order (right panel) Edgeworth and Gamma PDF expansions as functions 
of v = 8/cr for several values of <r, and Sj as labeled in the Figures. The dotted, dashed and continuous lines correspond to the Gaussian 
PDF, Gamma and Edgeworth expansions, respectively. 



value) for negative values of v = 5/a, and large modula- 
tions for v > 0. For some range of parameters and v > 
(ie. S > 0), the Gamma PDF also exhibits negative prob- 
abilities but they are typically much smaller than in the 
Edgeworth expansion. This problem was less apparent in 
Fig. w, when working at the first order of the series (left 
panel). 

The long-dashed vertical line corresponds to v = — 1/cr, 
which marks the range of positive density: p > or 
v > —1/cr. This line is crossed both by the Gaussian 
PDF and the Edgeworth PDF, indicating that these dis- 
tributions cannot be used as physical distributions above 
some small value of a. 

5. COMPARISON WITH SIMULATIONS 

As it is suggested by Eqs(p3),(B5h, differences in both 
expansions might be slight, specially around the peak of 
the PDF, as to first order both expansions are formally 
equivalent. Thus, which one best fits numerical results is 
a matter of careful analysis. The Edgeworth expansion has 
been shown to provide a very good approximation to model 
the PDF resulting from the weakly non-linear gravitational 
growth from small Gaussian initial fluctuations. We next 
compare both expansions with N-body simulations to see 
if we find significant deviations in their predictions. 

We measure the PDF in 10 realizations of SCDM sim- 
ulations, n = 1 and T = 0.5, with L = I8O/1- 1 Mpc and 
N = 64 3 particles and a s = 1 (Croft & Efstathiou 1994). 
We consider several smoothing radius which correspond to 
different values of the variance, a 2 . The errors correspond 
the rms standard deviation in the 10 realizations. 

Figure || corresponds to a ~ 0.2. As can be seen in the 
plot both expansions including the second order produce 
very similar results, within the error bars. Note also that 
they both are significantly different from the Gaussian re- 



sult (dotted line). Similar results are obtained for lower 
values of a. 




Fig. 2. — Comparison of gravitational simulations with the sec- 
ond order Edgeworth and Gamma PDF expansions as a function of 
v = 5/cr. We use as parameters the measured values of a, S3 and 
S4 as labeled in the Figure. The dotted, dashed and continuous 
lines correspond to the Gaussian, Gamma and Edgeworth distribu- 
tion expansions. The inset shows a detail around the peak in linear 
scale. 

When a is small, a < 0.3, using the measured N-body 
(non-linear) values of a and Sj as parameters for the PDF 
expansion, it yields very similar numbers, within the error 
bars, to the ones obtained using the corresponding lin- 
ear a and non-linear PT predictions for Sj. The latter is 




a, 







Fig. 3. — Deviations from the Gaussian PDF for both expansions and in N-body simulations (symbols). The lower panels displays results 
for the first order. The upper panels show the expansions including the third-order terms. The solid line is given by the Edgeworth series 
while the dashed one shows the Gamma expansion. The left and right panels show different ranges in v. 



in agreement with what was suggested by Juszkiewicz et 
al. (1995), Bernardeau & Koffman (1995). For larger val- 
ues of the rms fluctuation, a > 1, differences in the PDF 
when using perturbative or non-linear parameters become 
larger. 

In Figure we can see that for a ~ 0.3 the Gamma 
expansion seems to be in better agreement with the nu- 
merical results than the Edgeworth specially for negative 
values of v. It is seen that both expansions yield similar re- 
sults, except for the negative tail which is better recovered 
by the Gamma expansion again. The Gamma distribution 
seems to perform slightly better also around v ~ 1 — 5 as 
can be seen in Figure ra (right panel). For small v the 
Edgeworth series yields relatively large negative probabil- 
ities of the order P{y) ~ 0.01. Thus, in this case it would 
be better to use the Gamma expansion if our priority is 
a better overall behaviour. This needs not be rigorously 
generic and the situation could change if we explore a dif- 
ferent domain of the parameter space. It mainly depends 
on how large the values of Aj are for the case of inter- 
est. Also note that larger values of a require higher orders 
in the expansions, which could then change their relative 
performances. 

Figure illustrates what happens for larger values of 
a, where these expansions are expected to break down. 
Only the first order is shown, but the agreement does not 
improve for higher orders. In this regime, a ~ 1, the 
expansions do not work well, as expected, and one has 
to use some different approach to model the PDF. The 
short-dashed line shows the results of one such model: the 
non- linear Spherical Collapse model (from Fig. 7 in Gaz- 
tanaga & Croft 1999). Errobars in the simulations are 
comparable to the size of the symbols. The Spherical Col- 
lapse model for the evolution of density perturbations can 
be used as a local Lagrangian mapping to relate the ini- 
tial and evolved fluctuation. The only input required is 



the linear variance a and its slope 7 (eg see Fosalba & 
Gaztahaga 1998a, b, Gaztahaga & Fosalba 1998). Notice 
however that, while the Edgeworth and Gamma expan- 
sions could be used to model an arbitrary PDF, given its 
moments Sj, the Spherical Collapse is only intended to 
model gravitational dynamics of cosmic fluctuations. 




Fig. 4. — Comparison of Nbody results (symbols) with the 
Gaussian (dotted line), the Edgeworth series (continuous) and the 
Gamma expansion (short dashed line). Both expansions are given 
to first order. The long-dashed line following closely the simulation 
points corresponds to the spherical collapse model. 

6. DISCUSSION & CONCLUSION 

In this paper we have introduced the Gamma expan- 
sion, as an alternative to the well-known Edgeworth se- 



ries, to model the gravitational evolution of the density 
field on large-scales in terms of its lower order moments 
Sj = fcj//t2 (see Eq([|)). The search for an alternative 
expansion is motivated by the fact that the Gaussian PDF, 
which is used as the parent distribution for the Edgeworth 
series, is not strictly well-defined for describing positive 
variates, such as the density field. As a consequence, the 
expansion built out of the Gaussian exhibits undesirable 
properties, namely, it predicts negative probabilities and 
allows for negative densities. In §@ we have provided an 
independent derivation of the Edgeworth series based on 
the saddle point approximation to the Legendre transform 
of the PDF. 

In the case of the Gamma expansion (see §||) , the basis 
for the series is given by the generalised Laguerre polino- 
mials. These are the counterparts to the Hermite polino- 
mials which appear for expansions around the Gaussian, 
and allow for the construction of a consistent series in the 
perturbative limit (when the variance is small). The co- 
efficients in both expansions can be writen in terms of 

A,/ = (Sj — Sf)/J\, where the reduced cumulants of the 
PDF, Sj, appear as differences with the corresponding cu- 
mulans, Sf of the parent PDF (over which we are expand- 
ing), ie. S { f } = for a Gaussian PDF, or S^ = (J-l)! for 
the Gamma PDF (see Eq(§§) and Eq@). The Gamma 
expansion recovers all the contributions that appeared in 
the Edgeworth series, as functions of Aj, plus some cor- 
rective terms, for second order in a or higher. The con- 
vergence of the series is safe as long as the PDF develops 
a tail close enough to exponential, ie. when Aj are small. 
Gravitational clustering predicts the appearence of such 
exponential tails in the weakly non-linear evolution of the 
cosmic density field, at least for Gaussian initial condi- 
tions. Therefore, the expansion introduced in this paper 
should constitute a good candidate to properly model clus- 
tering on large scales. 

We have carried out a detailed comparison of the per- 
formance of the Edgeworth and the Gamma expansions 
in the perturbative regime with respect to cosmological 
N-body simulations with a < 0.4. We have found that 
they both yield a very similar agreement with numerical 
results around the peak of the distribution. The Gamma 
expansion provides a better general match to the PDF 
on the tails. In particular, the negative density tail mea- 
sured in N-body simulations is accurately recovered from 
the Gamma expansion, unlike the case of the Edgeworth 
series. Nevertheless, in general, the performance of the 
expansions depend strongly on the values of Aj. The 
agreement is better for the expansion which has its par- 



ent moments Sf closer to those of the PDF we want to 
model. In other words, the smaller the Aj, the better the 
behaviour of the expansion of the PDF. 

In summary, the Gamma expansion provides an inter- 
esting alternative to the Edgeworth series without intro- 
ducing any additional mathematical entanglement. Both 
expansions have the same inputs (the cumulants of the 
PDF we want to model) and outputs (the recovery of the 
full PDF) and similar expressions (see Eq(pl) and Eq(p5|)). 
The proposed Gamma expansion is better suited forde- 
scribing a real PDF, because it always yields positive den- 
sities and the PDF is effectively positive-definite. It is also 
possible to expent the Gamma expansion to a multivariate 
form, using a generalization of the Laguerre polinomials. 
Finally, we stress that many of the arguments presented 
here are rather generic and they might be useful when 
adressing the problem of modeling PDFs in other contexts. 
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